Topological minigap in quasi-one-dimensional spin-orbit-coupled semiconductor Majorana wires 
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The excitation gap above the Majorana fermion (MF) modes at the ends of ID topological superconducting 
(TS) semiconductor wires scales with the bulk quasiparticle gap E qp . This gap, also called minigap, facilitates 
experimental detection of the pristine TS state and MFs at experimentally accessible temperatures T <C E qp . 
Here we show that the linear scaling of minigap with E qp can fail in quasi-\D wires with multiple confinement 
bands when the applied Zeeman field is greater than or equal to about half of the confinement-induced bandgap. 
TS states in such wires have an approximate chiral symmetry supporting multiple near zero energy modes at 
each end leading to a minigap which can effectively vanish. We show that the problem of small minigap in 
such wires can be resolved by forcing the system to break the approximate chirality symmetry externally with a 
second Zeeman field. Although experimental signatures such as zero bias peak from the wire ends is suppressed 
by the second Zeeman field above a critical value, such a field is required in some important parameter regimes 
of quasi- ID wires to isolate the topological physics of end state MFs. We also discuss the crucial difference of 
our minigap calculations from the previously reported minigap results appropriate for idealized spinless p-wave 
superconductors and explain why the clustering of fermionic subgap states around the zero energy Majorana end 
state with increasing chemical potential seen in the latter system does not apply to the experimental TS states in 
spin-orbit coupled nanowires. 

PACS numbers: 03.67.Lx, 03.65.Vf, 71.10.Pm 



I. INTRODUCTION 

Ref. proposed a ID spinless p-wave superconductor 
as a platform for end-state Majorana fermions and topolog- 
ical quantum computation (see also Ref. |2|). It has been 
shown that a dimensionally reduced ID version of a recently 
proposed 2D Rashba-coupled semiconductor heterostructure^ 
can realize such a system for experimental investigations^]. 
The ID system has the advantage that the minimum excitation 
gap above the end-state MFs, the so-called minigap, scales 
with the induced quasiparticle gap E qp ~ 1 K in the nanowire. 
Such a large minigap (compared with A 2 /ep ~ 0.1/iK in 
2D TS systems^, for example, in Sr 2 Ru04, a potential chi- 
ral p-wave topological superconductor 7 ) allows the realization 
of the pristine TS state and MFs at experimentally accessible 
temperatures by applying a strong Zeeman field. It has been 
shown by explicit calculations 4 5 that the end-state MFs thus 
produced can be probed in local tunneling experiments that 
evince a zero-bias conductance peak. Ref. [ 8 ], which indepen- 
dently proposed the ID structure, and Ref. [ 9 ] proposed an ac- 
Josephson experiment to probe the nanowire MFs. An alter- 
native direct way of having a large minigap in a TS system is 
to use a 2D hybrid semiconductor- superconductor structure 3 
with a very small Fermi energy in the semiconductor so that 
the 2D minigap, i.e. A 2 /cf, is intrinsically large. 

The basic idea behind the dimensional reduction of the 
2D structure and the associated minigap is as follows 4 -! The 
2D (xy) -plane Rashba-coupled semiconductor in the TS state 
supports a single gapless chiral Majorana edge mode with dis- 
persion E(k y ) = (A/kp)k y on an edge parallel to, say, the y- 
axis. From the mathematical equivalence, H < j^ G (k y = 0) = 
Hg® G , where H^ G is the Bogoliubov-de Gennes (BdG) 



Hamiltonian of a wire along the x-axis, it immediately fol- 
lows that a ID nanowire along x must have a zero energy MF 
eigen- solution localized at the end (E(k y = 0) = 0). For 
minigap, note that the 2D Hamiltonian has no other eigen- 
solution on the edge other than the gapless Majorana mode 
itself. This implies that the dimensionally reduced ID Hamil- 
tonian should also have no other sub-gap solution other than 
the zero energy state at the end. Hence, the minigap in the ID 
problem should really be equal to the quasiparticle gap E qp 
induced in the nanowire. Viewed another way, in going from 
the 2D plane to the ID wire, as the width L y in the ^/-direction 
is reduced, the energies of the quantized edge modes scale as 
~ l/Ly, i.e. as the inverse of the confinement size of the wire 
in the transverse direction - we are assuming here that the 
size of the wire in the third z-direction is much smaller than 
that in the ^-direction, but what matters is simply the trans- 
verse confinement size being very small compared with the 
length of the wire. Since this quantity ultimately diverges in 
the strict ID limit L y — > 0, it implies the minigap to be equal 
to the quasiparticle gap E qp . The end- state MFs protected by 
such a large minigap can be probed in local tunneling experi- 
ments without having to worry about other low energy states 
with energies comparable to experimental temperatures. The 
zero energy modes should manifest themselves in zero-bias 
conductance peaks which exist only above a critical value of 
the Zeeman splitting required for the topological state. These 
conjectures, numerically confirmed for strict ID wires 4 , have 
appeared with calculational details in Ref. (5). 

In this paper we revisit the question of the minigap in ID 
wires and show that the linear scaling of minigap with E qp is 
valid only for strict ID wires (or wires where the transverse 
confinement induced band gap far exceeds all other energy 
scales but excluding the Fermi energy for multiband occu- 
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pancy). For quasi- ID or multi-band wires^ where the con- 
finement band gap is comparable to the applied Zeeman en- 
ergy seal J 11 * 12 * the above linear scaling can spectacularly break 
down. The parameter regime where the relevant confinement 
band gap is comparable to the Zeeman energy is the so-called 
'sweet- spot' regime^ 2 ^, designed to maximize the robustness 
of the topological state to inevitable spatial chemical poten- 
tial fluctuations arising, for example, from unintentional ran- 
dom background charged impurities in the semiconductor. We 
show that in this regime, because of an approximate chiral- 
ity symmetry^Sl f th e topological state, multiple near zero 
energy modes appear at the same end and the nanowire mini- 
gap effectively vanishes. The chirality symmetry in this case 
refers (see Sec. Ill) to the existence of a unitary operator that 
anti-commutes with the Hamiltonian. This symmetry, strictly 
valid for quasi- ID wires only in the absence of an interband 
Rashba coupling, ensures multiple zero energy MF modes 
on a given end of the Majorana wire. For a finite interband 
Rashba coupling the chirality symmetry is no longer exact, 
but we find that the quantitative effects of experimentally re- 
alistic interband Rashba couplings are insignificant and call 
the resultant Hamiltonian approximately chiral symmetric. 

Because of the possibility of multiple zero and near zero 
modes from the same end (with their total number given by 
the integer TV in Fig. [1]), we use here a slightly generalized 
definition of the minigap as explained below. If TV is allowed 
to take only values 0, 1 (i.e., if the quantum wire had no hid- 
den approximate chirality symmetry), the minigap would be 
strictly defined as the minimum excitation gap above the zero 
energy MF state. This definition of the minigap, however, is 
inadequate when TV is allowed to take any integer values in 
the relevant parameter space. The integer TV includes, with 
the hidden chiral symmetry of the nanowire weakly broken 
by the interband Rashba term, (TV mod 2) exact zero energy 
MF states and the remaining (TV-(TV mod 2)) near zero en- 
ergy states with very small energy splittings ~ 10~ 2 E qp , see 
Fig. [3]. This is because only in the presence of exact chirality 
symmetry an arbitrary number of MF modes can be supported 
by a given end of a Majorana wire. A finite interband Rashba 
coupling breaks the chirality symmetry, splits the energies of 
the zero modes in pairs, but does so only weakly. So in the 
presence of an experimentally realistic interband Rashba, we 
will show below that there are no exact zero modes for states 
with even TV, but there are TV near zero modes. Similarly, 
in states with odd TV, TV — 1 zero modes turn into near zero 
modes by the interband Rashba term, while exactly one mode 
remains a zero energy MF mode. The remaining single MF 
on a given end cannot further be split because it is protected 
by the particle-hole symmetry of the BdG Hamiltonian which 
cannot be broken in a superconductor. 1 

The slightly generalized definition of the minigap can be 
understood by noting that, when, say, the TV = 1 state 
(Fig. [1]) is being probed by a zero bias tunneling peak, for 
an unambiguous determination of the MF one must also en- 
sure that the tunneling in reality is not probing the near zero 
energy end-modes in the adjoining TV = 2 state. This implies 
that the definition of the minigap be expanded to include also 
the energy of the near zero modes in the parameter regime 



nearby to the MF. This minigap should also be raised suffi- 
ciently above the experimental temperatures to exclude un- 
intentional contributions from near zero energy states in the 
relevant parameter space. As we will show, the minigap, thus 
defined, effectively vanishes for the semiconductor nanowire 
because of the hidden approximate chiral symmetry. As in the 
case of 2D TS systems, it then becomes technically very diffi- 
cult (although may not be impossible 15 ) to probe the physics 
of isolated MFs and non-Abelian statistics at experimentally 
accessible temperatures^^! 

Fortunately, in ID, we find that there is a solution to the 
small minigap problem which involves forcefully breaking the 
nanowire chirality symmetry by applying a second Zeeman 
field transverse to the wire. The second Zeeman field, orthog- 
onal to the one realizing the topological state, enhances the 
energies of the near zero energy modes in pairs but leaves the 
MF state intact, before removing all TV states above a criti- 
cal value at which the system transitions into a topologically 
trivial state with TV = (Fig. [5]). Thus, if a zero bias conduc- 
tance peak from an end of the nanowire 5 splits with increasing 
the second Zeeman field, the splitting is given by the minigap 
created by the second Zeeman field. For the states with var- 
ious TV (see Fig. [1]), a zero bias conductance peak from the 
TV = 1 state (with no near-zero energy mode) should show 
no splitting with the transverse field. The zero bias peak from 
the TV = 2 state splits with an applied transverse Zeeman field 
because of the externally induced broken chirality symmetry. 
The zero bias peak in the TV = 3 state is expected to decrease 
in height with the applied field but (TV-(TV mod 2)=2) new low 
energy peaks are split off following the energies of the near- 
zero modes with the increasing second Zeeman field. Finally, 
all low bias peaks should disappear for sufficiently strong sec- 
ond Zeeman field as the system transitions into a topologically 
trivial state with TV = (Fig. [5]). 



II. HAMILTONIAN, PHASE DIAGRAM, AND MULTIPLE 
NEAR ZERO MODES IN QUASI-1D WIRES 

In a recent paper^Niu et ah have shown that a ID spinless 
p-wave topological superconductor can in principle support 
any integer (not just or 1) number of MFs at each end. This 
has prompted two of us^l to examine more closely the anal- 
ogous case of a ID Rashba-coupled semiconductor nanowire 
with Zeeman splitting proximity-coupled to a s-wave super- 
conductor. We have shown that this system - the so-called 
semiconductor "Majorana nanowire" - can also support any 
integer number of MFs at a given end under invariance of the 
chirality symmetry. Even though the parent 2D system 3 and 
others symmetry-related to H 6 I 21 I 22 I are in topological class D, 
we have clarified that the topological class of the ID semi- 
conductor under chirality symmetry is BDI. The class BDIis 
characterized by an integer Z bulk topological invarianfEEl 
and thus allows multiple MFs (with number equal to the in- 
teger invariant) even from a single end. That the topological 
class of the ID system should be BDI follows from the dimen- 
sional reduction arguments mentioned in the introduction^. If 
we think of the ID end MF of a nanowire as the dimensionally 



3 




V 



20 40 60 80 100 120 140 
Zeeman field T/E a 

FIG. 1. (Color online) Phase diagram of a quasi- ID nanowire with 
proximity-induced superconductivity as a function of the Zeeman 
field T and the chemical potential /i. Different superconducting 
phases are characterized by the number N of low-energy modes with 
E n « at each end of the wire. The energies of these modes identi- 
cally vanish, with N denoting the number of exact zero energy MFs 
at each end, for H nrn (a y — 0) (see Eq. ([TJ). For H nrn (ay — a), 
(N-(N mod 2)) modes are near zero modes with very small energy 
splittings rsj 10 -2 E qp where E qp is the induced quasiparticle gap 
~ 1 K. The near-zero energy modes are very robust to all reason- 
able perturbations including disorder except to a strong Zeeman field 
perpendicular to the length of the nanowire. With the representative 
parameters of this paper (E a — m*a 2 — 0.053 meV and an ef- 
fective g* ~ 50 for InApl), the N = 1 region starts at B x ~ 0.7 
T. Due to the existence of multiple robust near zero modes at each 
end of the quasi- ID wire the entire region beyond B x ~ 0.7 T is 
effectively gapless. 

reduced version of a 2D gapless edge MF mode, then since in 
2D there are in principle Z edge MF modes allowed (class D 
in 2D is characterized by a Z invariant) it follows that there 
must be Z end MFs allowed also in ID. Hence the nanowire 
ID system should also be with a Z, not Z2, invariant, which 
leads to BDI being the appropriate topological class. (A sim- 
ilar dimensional reduction argument should predict that the 
topological class of the Fu-Kane system 20 should reduce from 
Dili in 2D to D in ID both having a Z2 invariant.) To prop- 
erly define the BDI class and the associated Z invariant for 
the nanowire, we had to invoke 18 a hidden chirality symme- 
try of the ID system and show that, under invariance of this 
symmetry, the ID nanowire can support an arbitrary integer 
number TV (= Z) of zero energy MFs at each end. Below 
we make a departure from the case of strict ID nano wires and 
consider the experimentally realistic case of quasi-YD wires 
with Rashba and Zeeman couplings and proximity induced 
superconductivity. We show that it retains an approximate 
chirality symmetry and thus allows the realization of multiple 
near zero energy modes at each end in experimentally relevant 
parameter regimes. 

We consider a rectangular semiconductor (SM) nanowire 
with lengths L x ^> Ly ^> Lz proximity coupled at the in- 
terface z = to an s-wave superconductor (SC) with a su- 
perconducting gap A = 1 meV. The derivation of the low- 
energy effective Hamiltonian of the semiconductor wire that 



FIG. 2. (Color online) Confinement band configurations for dif- 
ferent values of the integer N. N gives the number of near zero 
energy modes at each end of quasi- ID nano wires with Hamilto- 
nian H nrn (a y — a) given in Eq. |TJ. (a) Band configuration for 
N = 0,1. n~ and n + denote a pair of sub-bands for each confine- 
ment band index n and move in opposite directions in response to a 
Zeeman field, (b) Sub-band configuration for higher Zeeman fields 
producing N = 2. Note that the sub-bands 1 + and 2~ are switched 
in energy by the Zeeman field, (c) Confinement band configuration 
capable of producing N = 3. (d) The experimental set-up consisting 
of a nanowire on top of a superconductor. The thick (red) arrows in- 
dicate applied magnetic fields. B x produces the TS state, while B y 
is an additional transverse field proposed here to externally break the 
chirality symmetry and produce a reasonable minigap. 

includes the induced s-wave pairing and the energy renormal- 
izations due to proximity to the SC is described in detail in 
Ref. 1 12|. For an infinite wire (infinite in the £ -direction), the 
BdG Hamiltonian of the multiband system has the form, 

H nrn (/c) [fnm(^) f^^nm] / T~z H~ ^^nm^x^z 

H~ Otkb ^nmO ~y7 ~z ^^yQum^x H - ^■nmO~yTy ) (1) 

where k = k x is the wave number in the x-direction, cr^ and 
T{ are Pauli matrices associated with the spin- 1/2 and particle- 
hole (p-h) degrees of freedom, respectively, and we have used 
the basis (u^,u^,v^,v^) for the p-h spinors. In Eq. (fTl n and 
m label different confinement induced sub-bands described 
by the transverse wave functions <j) n (if) oc sin(?27ri/ j Ly) , 6 nm 
describes the SM spectrum in the absence of spin-orbit cou- 
pling, [i is the chemical potential, and T = g* iibB x /2 is 
the external Zeeman coupling along the x-direction which is 
needed to create a topological state. In the presence of an in- 
homogeneous SM-SC coupling, the induced superconducting 
pairing A nm contains non- vanishing inter-band components. 
The inter-band Rashba coupling is described by the renormal- 
ized matrix elements q nm ex {4> n \d / 'dy\(j) m ) , which couple 
transverse states with opposite parity. Note that in the ef- 
fective Hamiltonian all energies are renormalized due to the 
proximity effect by a Z-factor that depends on the strength of 
the SC-SM coupling, as described in Ref. fV2\ . The effective 
parameters e nm , q nm , A nm are calculated numerically fol- 
lowing the procedure described in Ref. fT2ll . The value of the 
Rashba spin-orbit coupling used in the calculations is a = 0.1 
eV A. For the calculations of the phase diagram (see Fig. [1]), 
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the coefficient of the transverse Rashba coupling a y that pro- 
vides inter-band coupling is taken to be equal to a which is 
the experimentally relevant value. This term, however, breaks 
the exact chirality symmetry of H n m = H nrn (a y = 0) to 
only an approximate one for H nm . Thus, a finite a y turns the 
energies of some of the exact zero modes of H nm into only 
approximate zero modes for H nm . Despite this, as we show 
in detail below, the energy splittings of these (N-(N mod 2)) 
near zero modes due to a y = a are very small ~ 10~ 2 E qp . 
On the scale of E qp , therefore, the low-energy spectrum of 
H nm {a y = a) has the zero and near zero modes all show 
up as zero energy modes (see Fig. [3]), whose number (from 
each end) is given by the Z topological invariant calculated for 
Hnm in the appropriate parameter regime. Thus, in this paper, 
we will calculate the integer topological invariant (the number 
W, see Eq. ^) for H nm and compare it with the number of 
near zero modes from a given end calculated by diagonalizing 
the full Hamiltonian H nrn (a y = a). 

The Hamiltonian H nrn (a y = a) has been recently stud- 
ied extensively in connection with realizing MFs in quasi- ID 
nanowires 11 12 . Since the Hamiltonian explicitly takes into 
account multiple confinement band occupancy, it not only 
lifts the stringent condition of strict one-dimensionally of the 
nanowire but also allows higher carrier density in the topo- 
logical states. Even more importantly, for values of the Zee- 
man coupling about half the confinement induced band gap at 
k = 0, the system becomes topological over a wide range of 
the chemical potential /i, see Fig. [1]. The quasi-lD nanowire, 
thus, allows the widest possible variation of /i without cross- 
ing a topological phase transition, leading to the maximum 
robustness of the end-state MFs to spatial disorder. This can 
be most easily understood as follows: Each confinement in- 
duced band index n corresponds to a pair of spin-orbit bands 
which are degenerate at k — in the absence of T. A finite T 
lifts this degeneracy and the two sub-bands in a given pair, say 
{n + ,n - }, move up and down (in opposite directions along 
energy axis) in response to increasing T. When the chemical 
potential falls in the gap created between n + and n~ at k = 0, 
the system is topological and hosts a single zero energy MF 
at each end of the wire. However, with increasing values of 
r - \E n (k = 0,T = 0) - E n+1 (k = 0,r = 0)|/2 the en- 
ergies of the bands n + and (n + 1)~ come close, and it can 
be shown that 11 12 in this regime (the so-called 'sweet spot' 
regime) the system is topological over a much wider regime 
of ii than is allowed for smaller values of T. The wide vari- 
ation of fi allowed in the same topological state can be seen, 
for instance, from the TV = 1 region in Fig. [1], TV indicating 
the number of zero energy modes at each end. 

From Fig. [1] it is clear that, with increasing values of 
T at fixed fx (the solid horizontal line), the system makes a 
series of topological transitions at which the integer TV in- 
creases by unity. For a fixed T, a similar trend of discon- 
tinuous shifts of TV is also visible with increasing values of 
fx (the solid vertical line in Fig. [1]). Later we will explain 
the jumps of the integer TV in terms of discontinuous shifts of 
the Z topological invariant of the chiral symmetric Hamilto- 
nian Hnm- In the limit of vanishing proximity induced pair 



potential, the integer enhancements of TV with increasing T 
can be understood as follows: For values of T at which the se- 
ries of bands {n~,n+} are 'normal-ordered' (see Fig. [2a]), 
the only possibilities for TV (equal to the integer topological 
invariant W (eq. ^) of H nm ) are the integers and 1. For 
r > \Ei(k = 0,T = 0) -E 2 (k = 0,T = 0)|/2, the bands 
1 + and 2~ switch in the energy axis, see Fig. [2b]. In this 
configuration of the bands, for values of fi corresponding to 
two Fermi surfaces from the lowest two bands 2 _ ), the 
topological W number, hence TV, becomes equal to 2. The 
system thus can support two independent exact zero energy 
MFs at each end for a y = 0. As mentioned before, a fi- 
nite a y = a raises the energies of this pair of end states by 
a minute amount ~ 10~ 2 E qp , where E qp gives the energy 
of the next higher excited state in the nanowire. For the full 
Hamiltonian H n m{&y = ot), TV = 2 therefore indicates two 
near zero energy states localized at each end. Proceeding this 
way, with further increase of T, when the configuration of the 
subbands become similar to that in Fig. [2c] TV becomes equal 
to 3 .While TV = 3 indicates three exact MFs from each end 
for Hnm, for H n m{oiy = a) it indicates only one end-state 
MF and two other near zero energy localized states with ener- 
gies ~ 10~ 2 E qp . The appearance of multiple near zero energy 
states at each end of the nanowire effectively reduces the mini- 
gap to ~ 10 ~ 2 E qp . It is in this sense we claim that the linear 
scaling of the minigap with E qp , valid in the case of strict ID 
wires 4 5 or when the confinement band gap far exceeds all en- 
ergy scales excluding the Fermi energy, breaks down in the 
case of quasi- ID wires in the experimentally relevant sweet 
spot regimes. We will now explain the occurrence of multiple 
near zero modes at each end of the multiband wire in terms of 
an approximate chirality symmetry and a topological integer 
invariant hidden in Eq. ([if. 

III. CHIRALITY SYMMETRY AND Z INVARIANT 
DESCRIPTION OF THE PHASE DIAGRAM 

For a discussion of the chirality symmetry of the nanowire 
we first take a y = 0. We will show below that this transverse 
Rashba term that couples the confinement induced transverse 
bands indexed by n, m breaks the chirality symmetry and pro- 
duces a tiny minigap ~ 10~ 2 E qp . We consider the Hamil- 
tonian Hnm = H n m{&y = 0)- It can be seen by explicit 
construction that H n m anticommutes with a unitary operator 
S r x , 

{Hnm,S}=0. (2) 

Here, the 'chirality' symmetry operator S = r x can be writ- 
ten as the product of an artificial 'time reversal' operator JC 
and a particle-hole transformation operator A = r x • JC where 
JC is just the complex conjugation operator. It is easy to check 
explicitly that H n m commutes with the complex conjugation 
operator JC and anticommutes with the p-h transformation op- 
erator A, and hence it anticommutes with the chirality opera- 
tor S = JC • A = r x . The existence of all three symmetries - 
'time reversal', particle-hole, and chirality - ensures that H n m 
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FIG. 3. (Color online) Correspondence between the number of near- 
zero modes of H nm (ot y = a) and the Z topological invariant of the 
chiral symmetric Hamiltonian H nrn for a system with \i — 54.5 E a 
and different values of the Zeeman field T. Upper panel: Low- 
energy BdG spectra. The number of near zero energy modes at 
each end of the wire increases from N — (at V — lO^a) to 
N = 3 (r = 80.E a ). In addition, for certain values of the Zeeman 
field, there are localized in-gap states with T-dependent energies of 
the order of the bulk quasiparticle gap E qp . Lower panel: Depen- 
dence of 6(k) (Eq. ([4])) on the wave vector k. The angle 0(k) is 
defined as a continuous function of k with 0(ir) = 0. Note that for k 
much larger than the Fermi wave vector 0(k) is practically vanishes. 
Since 0(0) — 0(k) is an antisymmetric function of k on the interval 
— 7r < k < 7r, the winding number W is given by the angle at k = 0, 
W = 0(0) /tt. 



is in the BDI symmetry clas d 13 ^ 14 * characterized by an integer 
topological invariant which we call W . From Eq. ^ it fol- 
lows that the large square matrix Hamiltonian H nm can be 
off-diagonalized in a basis in which the unitary operator S is 
diagonal: 



UH nm (k)U^ 



A(k) 
A T (-k) 



(3) 



Following Ref . EH we now define the variable, 

z(k) = exp(»0(fc)) = Det(A(k))/\Det(A(k))\, (4) 
and calculate the integer invariant, 

rk= * dz(k) 



W 



zi [ 

f Jk=o 



z(k) 



(5) 



which is an integer (W G Z) including zero. 

Note that in the presence of a finite transverse Rashba cou- 
pling, H nrn (a y ^ 0) does not anti-commute with r x . Hence, 
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FIG. 4. (Color online) Dependence of the quasiparticle gap and the 
integer topological invariant W on the Zeeman field (upper panel) 
and the chemical potential (lower panel). The parameters correspond 
to the horizontal and vertical cuts through the phase diagram indi- 
cated in Fig. [1]. The winding number W is calculated in the absence 
of transverse Rashba coupling (a y = 0), while the quasiparticle gap 
is determined using the full Hamiltonian H nm (ay — a) in Eq. ([TJ. 
Note that the phase boundaries can be identified by either the van- 
ishing of the gap or the discontinuities of the winding number. The 
inset (upper panel) shows the behavior of 0(k) near the topological 
phase transition at \i — 54.5^ a and V « 21. 55Ea. Note that along 
the phase boundaries 0(k) becomes ill-defined at k = 0. 

the Hamiltonian matrix is no longer off-diagonalizable in the 
diagonal basis of S and the number W cannot even be defined. 
A finite a y thus breaks the chirality symmetry. However, even 
though the invariant W cannot be defined for a finite a y , since 
the experimental value of a y = a makes only a minute contri- 
bution ~ 10 -2 E qp to the energies of the near zero energy end 
states (compared to the energies ~ E qp of the next higher or- 
der excitations), we calculate the Z invariant for B. nrn (i.e., by 
taking a y = 0) and compare it with the phase diagram derived 
by computing the low energy spectrum of the full Hamiltonian 
H n m(oty = ol) in Eq. |T|). The results are shown in Fig. [3]. It 
is clear that the different topological phases of the full Hamil- 
tonian H nm (a y = a) characterized by different numbers TV 
of near zero energy end states can be characterized by differ- 
ent values of the integer W calculated for the corresponding 
reduced Hamiltonian H nm (a y = 0). 

The accurate correspondence of the topological quantum 
phase transitions separating phases with different values of 
N with the change of the integer invariant W is shown even 
more clearly in Fig. [4]. In this figure, we have plotted the in- 
duced quasiparticle gap E qp with the tuning parameters T (up- 
per panel) and the chemical potential fi (lower panel). Phases 
with different values of N are separated by points in the pa- 
rameter space where E qp vanishes. As is clear from Fig. [4] 
these are also the points at which the integer W changes dis- 
continuously, and the value of W for a finite value of E qp 
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FIG. 5. (Color online) Upper panel: Dependence of the minigap on 
the transverse Rashba coupling a y for \i — 54.5E a and two different 
values of the Zeeman field in the phase with TV = 2. Note that 
the low-energy modes have a finite energy E\ — E 2 that increases 
approximately linearly with the transverse Rashba coupling strength. 
For a y = a the minigap attains values ~ 1 — 2 x 10~ 2 E qp (see 
Fig. [4] for comparison). Lower panel: Dependence of the lowest 
energy states on a transverse Zeeman field for \i — 54.5E a and T = 
bOEa. Note that E± = E 2 if a y = 0, while E\ / E 2 in the presence 
of both the transverse Zeeman field T y and the transverse Rashba 
coupling dy. For T y > 20 E a the quasiparticle gap itself collapses 
(black line) and the nanowire makes a transition to a topologically 
trivial phase with no MFs. 

corresponds to the number TV of the near zero energy modes 
localized at each end of the quasi- ID nanowire. 



IV. CHIRALITY BREAKING AND NON-ZERO MINIGAP 
WITH TRANSVERSE ZEEMAN COUPLING 

As we have seen above, despite breaking the chirality sym- 
metry, the quantitative effects of a y = a on the energies of the 
localized end modes of quasi- ID nanowires are minimal. In 
particular, even though the state with TV = 2 no longer has an 
exact zero energy mode localized at the ends, the energy split- 
ting caused by a y is very small, ~ 10~ 2 E qp (E qp ~ 1 K). The 
situation is similar for TV = 3 where a y = a raises the ener- 
gies of two of the zero modes to ~ 10~ 2 E qp while the third 
one remains at exact zero energy which is a Majorana mode. 
We emphasize that the lowest energy state for TV = 3 is a Ma- 
jorana mode. Any finite gap acquired by this state is entirely 
due to finite size effects. 24 The finite transverse Rashba cou- 
pling thus creates a small minigap above the end state MFs, 
but such a small gap is difficult to resolve experimentally at 
experimentally realistic temperatures. 

We propose to solve the small minigap problem of quasi- ID 
nanowires by applying an additional transverse Zeeman field 
T y = g* [isBy/2 in addition to the longitudinal one needed 
to create the TS state itself (see Fig. [2d]). Note that such a 
Zeeman field is still parallel to the plane of the nearby bulk 
superconductor and thus does not create additional problems 



with the s-wave superconducting pairing. In the presence of 
this term, the BdG Hamiltonian of the nanowire becomes, 

H nm( k ) = H nm(k) + T y 5 nm <T y . (6) 

It can be easily checked that both terms with coupling con- 
stants T y and a y retain diagonal elements even in the basis in 
which S = r x is diagonal. This implies that H' nrn is no longer 
chiral symmetric with S as the chirality operator. Based on 
the available discrete symmetries (time reversal, particle-hole, 
complex conjugation, etc.) it is not possible to find any unitary 
symmetry operator that anticommutes with H' nrn . It follows 
that the chiral symmetry hidden in H nm has now been bro- 
ken externally by T y . In the absence of T y if the number of 
(near) zero modes is even, the transverse field creates a gap 
for all of them resulting in no MF edge mode. In cases where 
the number of zero modes is odd for Y y = 0, the transverse 
field opens a gap for all of them save one, and the system is 
reduced to having only one non-degenerate MF edge mode (at 
each end). As a consequence, the minigap above the end-state 
MFs is now externally tunable by tuning T y . It should be kept 
in mind, however, that T y cannot be increased to arbitrarily 
large values to attain a large minigap ~ E qp . This is because, 
above a critical value T y = Ty, the quasiparticle gap E qp 
itself closes and the system makes a transition into a topolog- 
ically trivial state with no MFs (see Fig. [5]). Keeping these 
constraints in mind, we show in Fig. [5] that a reasonable T y 
(< Ty) can raise the minigap to about ~ 0.1E qp ~ 100 mK. 

For experimental signatures such as the zero bias conduc- 
tance peak from wire endP, our analysis predicts that in the 
sweet spot regimes for T y = the zero bias peak in the TV = 1 
state of Fig. [1] is expected to continue also in the TV = 2 
state because of the existence of two near zero modes from 
the same end. In the presence of increasing T y , the peak in 
the TV = 1 state will disappear above T y > Ty. The zero 
bias peak in the TV = 2 state, however, is expected to first 
split into two with increasing T y ( because of the splitting 
(Fig. [5]) induced by the broken chirality symmetry), before 
disappearing altogether for T y > Ty in the topologically triv- 
ial state. Thus, just for the observation of the zero bias peak 
from the wire ends, a non-zero T y is not required. However, 
to unambiguously probe the physics of isolated MFs and non- 
Abelian statistics in the sweet spot regimes, an optimum trans- 
verse Zeeman field is required to raise the minigap sufficiently 
above the experimentally accessible temperatures. 

V CONCLUSION 

The minigap above the end- state Majorana fermions in ID 
topological superconducting nanowires scales linearly with 
the bulk quasiparticle gap E q p . Thi s fact, now confirmed 
in many numerical calculation j 4 l 5 l 12 l 23 J solves the problem of 
small minigap (~ A 2 /e^ ~ 0.1/iK) of 2D TS systems where 
MFs appear in bulk order parameter defects such as vortices. 
Here we show that the linear scaling of minigap with E qp is 
valid only for strict ID wires or where the confinement in- 
duced band gap far exceeds all other energy scales exclud- 
ing the Fermi energy. We show that the linear scaling fails 
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in parameter regimes (called 'sweet spot' regime) where the 
confinement band gap is in the range of the Zeeman energy 
scale for which the MFs have been proposed 11 12 to be most 
stable against spatial disorder and chemical potential fluctu- 
ations. The failure of the linear scaling arises from an ap- 
proximate chiral symmetry of quasi- ID nanowires leading to 
the possibility of multiple near zero energy modes from each 
end. The near zero modes are robust to all reasonable per- 
turbations including disorder, except to an additional trans- 
verse Zeeman field that can be used to break the approximate 
chiral symmetry externally. We show that the nearly vanish- 
ing minigap of quasi- ID nanowires in this regime can be re- 
stored to experimentally accessible temperatures by applying 
the second Zeeman field perpendicular to the length of the 
nanowire. For experimental signatures such as the zero bias 
conductance peak from the wire endP, the additional Zeeman 
field removes the zero bias anomaly in the topologically triv- 
ial phase for T y > Ty. An optimum T y is needed, however, 
to probe the topological physics of isolated MFs in sweet spot 
regimes because it can raise the minigap above the experimen- 
tally realistic temperatures. 

The present paper provides the theoretical description of 
the minigap for the experimentally relevant quasi- ID semi- 
conductor quantum wires with strong spin-orbit coupling in 
proximity to a superconductor and external Zeeman fields. 
Previously reported calculations for the minigap 25 26 based 
on idealized quasi- ID spinless p-wave superconducting wires 
do not apply to the experimental semiconductor nanowire 
systems where the interplay of spin-orbit coupling and Zee- 
man splitting leads to a helical topological superconductor in 
the presence of the proximity effect. Even though the two 
systems have some superficial similarities, notably the topo- 
logical class which is BDI under chirality symmetry^, the 
calculation of the minigap is one example where the differ- 
ence is explicit as explained below. For a quasi- ID spin- 
less p-wave superconductor Ref. |25 ] finds a large number of 
end-localized subgap states with very low energies clustering 
around the zero energy end state when the chemical potential 
crosses many confinement induced subbands. The interpreta- 
tion of this effect is in terms of each subband contributing its 
own MF state, which are then weakly split by the transverse 
couplings leaving only one exact zero energy eigenvalue when 
the number of relevant subbands is odd. We want to empha- 
size that this effect, which drastically reduces the minigap in 
spinless p-wave wires with many occupied subbands, is com- 



pletely absent in the spin-orbit coupled semiconductor wires 
at low applied fields. For the latter system it is clear from 
Fig. [3] (upper panel, T — 30^) that despite many subband 
crossings the energy of the first excited state above the MF 
end state for N = 1 is close to the bulk quasiparticle gap E qp . 
Thus, in this case, there is no effect akin to the crowding near 
zero energy effect discussed in the context of spinless p-wave 
superconductors. Therefore, the fermionic near-zero energy 
states of Ref. l25l are not a problem for the topological phase 
in real semiconductor nanowires although they are generically 
present in the Kitaev-type idealized model of spinless ID p- 
wave superconductor. Thus, the experimental semiconductor 
nanowires are in fact better systems for realizing Majorana 
modes than the idealized spinless p-wave ID system often 
studied in theoretical analysis. 

This crucial difference between the spinless p-wave case 
and the case of a semiconductor quantum wire with spin-orbit 
coupling can be understood with the help of the schematic 
energy level diagrams in Fig. [2]. Note that the effective 
p-wave superconducting order parameter in the basis of the 
semiconductor bands comes with a sign difference between 
the two bands in a given spin-orbit pair {n~, n + }. This ef- 
fect underlies the fact that with increasing /i, but staying at a 
fixed T, no extra near zero modes are generated in the wire 
ends despite more and more subbands being populated with 
increasing chemical potential. This is no longer true if the 
Zeeman coupling is large enough to switch the order of the 
subbands (see Fig. [2]) and only bands with one sign of the 
effective superconducting order parameter are populated with 
increasing /i. In this case, additional MFs will appear with 
increasing number of the occupied subbands, which will then 
be weakly split by the chirality breaking fields we discuss in 
the paper. While this is the case with the spin-orbit coupled 
quantum wires, the quasi- ID spinless p-wave superconductor 
has the same sign of the order parameter in all its confinement 
bands. This places the ideal ID p-wave system considered in 
Ref. [25 ] in the T = oo limit of the spin-orbit coupled quan- 
tum wires and only in this unphysical limit the minigap is re- 
duced in the experimental quantum wires by the multiband 
effects discussed in the context of spinless p-wave supercon- 
ductors. This problem obviously will not arise in the N = 1 
state (see Fig. [1]) of the real experimental systems which are 
spin-orbit coupled semiconductor nanowires. 

This work is supported by DARPA-MTO, NSF, DARPA- 
QuEST, JQI-NSF-PFC, Harvard Quantum Optics Center, and 
Microsoft-Q. 
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